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Abstract. This paper studies the scattering matrix §(S; h) of the problem 

-h 2 ip"(x) + V(x)tp(x) = Eip(x) 
for positive potentials V G C°°(R) with inverse square behavior as x — > ±00. It 
is shown that each entry takes the form E>ij(E\ h) = §ry(E; K)(l + fkTij(E; h)) 
where (E; K) is the WKB approximation relative to the modified potential 

V(x) + ^-(x)~ 2 and the correction terms Oij satisfy \dgO~ij (E; K)\ < C^E~ k 
for all k > and uniformly in (E, K) G (0, -Eo) X (0, Ho) where Eo, ho are small 
constants. This asymptotic behavior is not universal: if —h?d^ + V has a 
zero energy resonance, then S(E; H) exhibits different asymptotic behavior as 
E — » 0. The resonant case is excluded here due to V > 0. 

1. Introduction 

This paper revisits the much studied problem of determining the reflection and 
transmission coefficients for semi-classical operators of the form 

(f.l) P(x 1 hD):=-h 2 ^ + V(x) 

where V is real-valued and assumed to decay at infinity. There are two atypical 
features of this work, at least relative to the existing literature on this topic: 

(i) we wish to understand the zero energy limit, in fact uniformljQ in small h 

(ii) the smooth potential V decays like an inverse square at both endfl 

We remark that (i) and (ii) are closely related. Indeed, the (x}~ 2 decay is "critical" 
with respect to the zero energy limit in the sense that (x)~ 2 ~ s is easier and behaves 
very differently. In the semi-classical literature it is more customary to encounter 
the criticality of the Coulomb decay {x) ; the reason for this is that the Coulomb 
decay is critical for positive energies. Note that the numerology around these decay 
rates applies to all dimensions and not just to one dimension. The motivation 



The first author was supported by the National Science Foundation DMS— 0406193, DMS— 
0601226, DMS-0600369, the second by DMS-0653841, and the fourth by DMS-0405837. The 
second author thanks Fritz Gesztesy and Gerald Teschl for their interest in this work and comments 
on a preliminary version. 

^More precisely, the asymptotic analysis is carried out up to multiplicative errors of the form 
1 + O(h) where the O(h) needs to be uniform in small energies. 

2 The methods of this paper also apply to the case where the potential exhibits inverse square 
decay as x — > 00 and some other decay as x — > —00; for that, one of course needs to be able to 
carry out the scattering theory on x < 0. If the decay is |x|~ a with < a < 2, then 25 applies, 
whereas for a > 2 one can use classical scattering methods. 
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for considering this particular problem comes from several sources. First, smooth 
potentials which behave as an inverse square at one or both ends arise in several 
contexts in physics and geometry, for example in general relativity in connection 
with Schwarzschild and de-Sitter spaces, see [6]. Second, this paper is part of the 
program initiated in [23) and [24] . In fact, the analysis carried out here is an 
essential part in the solution of the "large angular momentum" problem from [24] . 

Let us briefly review some elementary features of scattering, cf. [7] and |17j : For 
simplicity, let h = 1 for now and write H = P(x, D). Recall that the Jost solutions 
f±(x; A) are required to satisfy 

Hf±(;X) =A 2 /±(-,A), f±(x,X)~e ±iXx as x ±oo 

Provided V € L 1 and A 7^ they exist and are uniquely determined as solutions of 
the Volterra equation 

(1-2) f + (x,X)=e^ + £° SHXi l~ X)) V(y)U(y,X)dy 

and similarly for /_(•, A). The resolvent kernel of H can now be expressed in the 
form 

iv c\2 , -nw-it \ f+(x,X)f-(y,X) f + (y,X)f-[x,X) 

(if - (A + i0)) (x, y) = X [ X >y] + X[x<y] 

for all A > where W(X) = W (/+(•, A), /-(■, A)). The reflection and transmission 
coefficients are defined by the relations 

MA)J+(, A) = r+(A)/_(-, A) + /_(•, -A) 
t_(A)/_(-,A)=r_(A)/ + (-,A)+/ + (-,-A) 

and satisfy 

(1.3) t-=t+, 1 = |<+| 2 + |r+| 2 = |L| 2 + |r_| 2 , r_=-f+*/* 

For fixed A > 0, consider the following bases of the space of solutions to the equation 
Hf = A 2 /: 

(/+(., A),/_(-, A)), (/+(-,-A),/_(-,-A)) 

The former is referred to as outgoing and the latter as incoming. In that case the 
matrix §(A) which transforms the coefficients of a solution relative to these bases 
satisfies 

't(X) r-(X) 
r+(A) t(X) 

It is called the scattering matrix and is unitary. Of special interest to us is the 
behavior as A — > 0+. Note that if 



S(A) = 



(x)\V(x)\ dx < 00 

then f + (x,X) — > f+(x,0) as A — > where the latter satisfies the limit equation 
of IlOI). viz. 



/+(a:,0) = l+ / (y-x)V(y)f + (y,0)dy 

J X 

It is known that S(A) is continuous in A > under this moment condition, see [12] . 
To describe the possible values of S(0), recall that H has a zero energy resonance 
iff /±(-,0) are linearly dependent or, equivalently, iff W(0) = 0. Furthermore, since 
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t(X) = — p^jy this is equivalent to t(0) 7^ 0. In conclusion, if zero energy is not 
resonant, then 



whereas in the resonant case 



S(0) 



5(0) 



-1 
-1 



t 



for some real r,t 6 [-1,1], t^O. 

If (x)V(x) £ i 1 (M), then the behavior of S(A) as A — > is completely different. In 
this paper, we focus on the border line case of positive inverse square potentials 
for (jl.ljl and H small (for the remainder of he paper, we now let H be a small positive 
quantity). It is precisely this case which arises in the geometric problem considered 
in [23], [23]. Our main theorem is as follows. We denote the energy by E = A 2 > 0, 
see above, and the scattering matrix of (jl.ip by 



i{E-n) = 



t(E;h) r-(E;h) 
r+(E;h) t(E-K) 



§2i(E;h) S 2 2(E;H) 



In view of (|1.3| it suffices to describe the first row of this matrix. In this paper, 
O(-) terms will be differentiablc functions and we will typically state bounds on 
their derivatives with regard to the relevant variables depending on the context. 

Theorem 1. Let V £ C°°(K) with V > and V{x) = [i 2 ± x~ 2 + 0(x" 3 ) as x -> ±00 
where [i+ ^ 0, /Lt_ ^ and d%0(x~ 3 ) = 0(x~ 3 ~ k ) for all k>0. Denote 



(1.4) 



V (x;h) : =V(x) + -(x) 



and let Eq > be such that for all < E < Eq and < K < 1, Vq{x] K) — E has a 
unique pair of solutions, which we denote by X2(E; K) < < X\(E] K). Define 

r x 1 {E;h) 

Edy 



(1.5) 



S(E;K) := / y/V (y;h) 

Jx 2 (E-h) 

poo 

T + (E;h):= Xl (E;h)VE- (V E - V (y; h) - Ve) dy 

Jx 1 (E-h) 



T-(E;h) := -x 2 (E;H)Ve - 



(E-h) 
x 2 (E-h) 



(^E-V (y;h)-VE)dy 



(1.6) 



as well as T(E; H) := T + (E; K) + T-(E; K). Then for all < h < Ho where ho = 
ho(V) > is small and < E < E 

S U (E; H) = e -mm+iT(E;K)) {1 + h au {E . h)) 

S 12 (E;H) = -ie-Tt T +( E W(l + ha 12 (E;h)) 
where the correction terms satisfy the bounds 

(1.7) \d k E o- 11 (E;h)\ + \d k E a 12 (E;h)\<C k E- k V k > 0, 

with a constant Ck that only depends on k and V . The same conclusion holds 
if instead of (fP|) we were to define V as V := V + h 2 Vi with V\ € C°°(R), 
Vi(x; h) = \{x)- 2 + 0{x- 3 ) as x -► ±00 with d k O{x- 3 ) = 0{x- 3 - k ) for all k > 
and uniformly m0<fi«l. 
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The addition of \{x)- 2 to V(x) is crucial and similar to the "Langer modifica- 
tion", see for example [§]; indeed, if we were to use V instead of Vb in Q1.5p . then 
the bounds (| 1 . T[) would fail due to a factor of log-E as E — ► 0. This is in contrast 
to potentials decaying like |a;| _a with < a < 2 for which the modification is not 
needed, i.e., the usual WKB ansatz works, see [25]. On the other hand, note that 
as long as Eo > E > e > the turning points Xj (E; ft) will remain bounded and 
the distinction between Vb and V is therefore moot. Indeed, the effect of passing 
from V to Vb and vice versa is merely a harmless factor of the form 1 + 0(h) where 
the O(-) term of course depends on e. In the range Eq > E > e > Theorem [T] is 
well-known and classical. See for example Chapter 13 of [19] as well as Ramond's 
work |22j for a more recent reference (Ramond, however, is more concerned with 
the scattering problem for energies close to the maximum of a barrier and he also 
assumes that the potential is dilation analytic). 

We remark that the infinite differentiability assumption on V can be relaxed to 
some finite amount of smoothness (in which case we can only ask for correspondingly 
many derivatives with respect to E), but we do not elaborate on this issue here. 
A more substantial problem is that of relaxing the positivity assumption. We 
conjecture that V > can be replaced by the strictly weaker assumption that 
zero energy is not a resonance of P(x,hD). Recall the definition of a zero energy 
resonance in this context, cf. [3], [25], and Section 3 of [24]: it means that the two 
subordinate zero-energy solutions at ±oo are linearly dependent (a "subordinate 
solution" at either end refers to the nonzero solution of P(x,hD)f = with the 
slowest possible growth at that end; it is unique up to a nonzero scalar factor). 

Note, however, that some condition is needed in Theorem [1] indeed, in [24] it 
was shown that for operators of the form considered in Theorem [1] with fjq_ = fjr_ = 
v 2 — j , v > | , and ft = 1 

(1.8) W(E;h) ~ E?-"(Wo + 0(E £ )) as E -> 0+ 

for some Wo ^ and e > provided there is no zero energy resonance. In the 
resonant case, it was shown in [21] that Wo = 0. The following relation between 
§n in ||OJ) and the Wronskian W(E; K) 

allows one to deduce (|1.8|) with Wo ^ from Theorem[T](note that for inverse square 
potentials S(E;h) behaves like |logi?| so that the apparent exponential behavior 
in (|1.6p turns into a power-law in E). This deduction also proves that Theorem [1] 
necessarily fails in the presence of a zero energy resonance. Another aspect of (|1.8p 
concerns the case of large h, say ft = 1. Indeed, it shows that Theorem [T] gives the 
correct behavior of the scattering matrix even in that case, but then the energy 
takes over as the small parameter. 

This paper is organized as follows. Section [5] constructs a fundamental system 
of zero energy solutions to (jl.ip via the usual WKB ansatz but for Vb rather than 
for V. Since we require uniform bounds in Theorem [T] as E — > 0, the construction 
of Jost solutions for positive energies which is carried out in Section [3] needs to yield 
the zero energy solutions in the limit E — > 0. We choose to reverse this process 
and show that Vb is precisely the right potential to use in the WKB method at zero 
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energy. The logic is simple: the WKB ansatz 



satisfies an equation of the form 



(-h 2 d 2 + v)^± = h 2 {- x - 2 /A + o(x- 3 ))^o,± 



where the x _2 /4 term on the right-hand side is universal for all potentials that have 
an inverse square decay as x — > oo as specified in Theorem [TJ Since this term has 
the same decay as V we need to bring it to the left-hand side leading to our choice 



The main technical work of this paper is carried out in Section [3] It is here 
that the (semi-classical) Jost solutions are constructed for all energies in the range 
< E < Eq. We use Langer's method which is based on the Liouville-Green 
transform, see Chapters 6 and 11 in [19]: switching to the new independent variable 



and to the new dependent variable w(Q = VC 7 / reduces P(x, hD)f = Ef, see (jl.ll) . 

to an Airy equation perturbed by a potential of size ti 2 . It is here that V > 
becomes relevant: it ensures that for all small E > there is a unique turning point 
x\(E) > and that Vo(x;h) > E for all < x < xi(E;K). Hence we can cover 
x > by the intervals £(0, E; H) < £ < and £ > 0. In each of these intervals we 
solve the perturbed Airy equations up to multiplicative errors of the form 1 + 0(h) 
where the O(-) term is uniform in E. It is in the range £(0, E; K) < £ < that the 
choice of Vq (rather than V) becomes decisive; this of course is to be expected as 
this range turns into the whole interval x > as E — » and WKB applied to V 
instead of Vq fails at E = 0, see Section [2] Theorem Q] is proved in Section [5] by 
evaluating the Wronskians 



at x = 0. Section [5] discusses the range of validity of Theorem [T] as the energy 
increases towards a unique non-degenerate maximum of a barrier potential. Fi- 
nally, the appendix describes a certain "normal-form" reduction of (jl.l| to a Bessel 
equation on a region containing the turning point. Even though we do not base 
our asymptotic analysis on this reduction (but rather the Airy equation), we still 
believe that this is of independent interest. 

Needless to say, there is a vast literature related to the semi-classical analysis of 
the Schrodinger equation and it is impossible to do any justice to it here. Somewhat 
curiously, however, there does not seem to be any literature on potentials which 
are globally smooth on the line and which exhibit inverse square decay. On the 
other hand, potentials which are exactly inverse square are of course ubiquitous, 
especially in the physics literature. For a recent paper in this direction involving 
WKB see [9] and for a time-dependent analysis see the recent papers [20] , [21] , as 
well as [4], [5] and the references cited there. Potentials which decay of the form 
\x\~ a , < a < 2, have been studied with similar objectives as here, see [13], [T8] 
and [25]. For other work on low energies see [2], [3], [8], and [26], as well as [ID] . 



oiV . 




2 



x>0 




W(U(;E),f_(;E)) 



W(f + (;E),f_(.,E)) 



6 OVIDIU COSTIN, WILHELM SCHLAG, WOLFGANG STAUBACH, SALEH TANVEER 



2. Zero energy solutions 

In order to motivate the choice of Vb in Theorem [I] we will now obtain a funda- 
mental system for the equation 

(2.1) -h 2 f'(x)+V(x)f(x) = 

on the half axis x > xq. Here we assume that V(x) = fi 2 x~ 2 + 0(x~ 3 ) with /i > 
as x — > oo and xq is chosen so large that V(x) > for x > Xq. As before, we require 
d x O{x~ 3 ) — 0(x~ 3 ~ k ) for all k > 0. Ignoring the O(-) term, we have on the one 
hand 

(-h 2 d 2 x + f i 2 x- 2 )xi ±a = 0, a 2 = j+n 2 h- 2 
On the other hand, with Q(x; H) := p?x~ 2 + h 2 x~ 2 /4, 



with some c ^ 0. This motivates the following result. 

Proposition 2. On x > xq a fundamental system of solutions for (12. ip is given 
by 

(2.2) ipj(x;h)=ipj(x;h)(l + ha,j(x;H)), j = 1,2 

with 

$i(x;H) =V (x;ti)^e* s ( x <V 
Mx;h)=V (x;h)^e^ s ^ h \ 



where V (x; K) = V{x) + |(i)- 2 , S(x; h) = J* q y/V (t;K)dt and 

(2.4) sup Id^ajfah)] < Ci^x' 1 

o<h<i 

for x > Xq, j — 1,2 and £ = 0, 1. Their Wronskian satisfies 

(2.5) W ^ 1 ^ 2) = zl {l + (h)). 
as H — > 0. 

Proof. Let us consider the case of i/'i- Hence, we need to find a\ so that "01 is a 
solution to the differential equation 

(2.6) - h 2 u" (x) + V{x)u{x) = 0. 

Substituting the first expression of (12.21) into the differential equation (12. 6|) yields 

(2.7) -h 2 [4>'( (1 + ftai) + 2nV4«i + T$ia"] + V^i(l + fioi) = 

Setting V2 := 5(2^ 2 - + and observing that -h 2 ^'{ + Vifi = -h 2 V 2 4>i, 

we deduce after dividing the equation by ipi 

(2.8) -(l + ha 1 )V 2 =h(a'( + 2-y-a' 1 ), 

0i 

We now note the following essential feature of V 2 (which was the reason for defining 
Vo as above): 

\V 2 (x)\ < Cx~ 3 , \d k x V 2 {x)\ < C k x- 3 - k Vfc>0 



ZERO ENERGY SCATTERING 7 

To solve (|2.8j) we multiply both sides by ip 2 and obtain 

(2.9) {a'^l)' =^-V 2 ${-a 1 V^l. 
Integration and using the definition of the ip± yield 

a[(x) = -j V 2 (y)^ 2 (x)^j(y) dy + J ai {y)V 2 (y)^ 2 {x)^l{y) dy 

(2.10) =~J V (x)iv (y)^-e^ s ^- s ^V 2 (y)dy 

Xq 

V (x) 5 Vb(y) ^e^ s ^- s ^V 2 (y) ai (y) dy. 

Strictly speaking, a\ — a\{x, h) but we suppress the h from the notation here. After 
integration in (|2. 10[) we obtain 

ax(x)=— I I Vo(x , ^Vo(y)-e-^ s ^ s(x ^V 2 (y)dydx' 



X J X 



(2.H) 



x J X 1 

%a ry 



-1 



x J x 



V (x')iv (y)^-e^ s ^- s ^V 2 (y)a 1 (y)dydx' 

V (x')iVo(y)^-e^ s W- s Wv 2 (y)dx'dy 
Vo(x'^V {y)^e^ s ^- s ^V 2 (y)a 1 (y)dx'dy. 



Furthermore, 



(2.12) f V (x')h^ s ^dx' = - \[e^ s ^ - e^ s ^}. 
J x 2 

and 

(2.13) V (y)^e^ f" V (x')h^ s ^dx' = —V (y)^ [1 - e f(S(v)-S(*))]. 

ii 2 

From this it follows that 

(2.14) ai{x) = \ f V Q {y)^[e^ s ^- s ^ -l]V 2 {y)[l + h ai {y)]dy. 

This is a standard Volterra equation. To solve it, we first introduce a new function 
p(x) given by 



(2-15) p{x) = / \V (y)-V 2 (y)\dy. 

Jxq 

In view of the decay of V 2 we see that the integrand here decays like y~ 2 so that 
p G L°°(xq, oo). Then we define a sequence etf (x), s — 0, 1, . . . , with a1(x) = 
and 

(2.16) at(x) = \ f X V (y)^[e^ S(y) ~ S{x)) ~ \]V 2 {y)[l + ha^ 1 {y)] dy 



We claim that 

(2.17) \at(x)~ar\x)\ - 
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To prove this we proceed by induction and observe that 

S(y)-S(x) = - ( V h (t,h)dt 



and hence for < x < y < x, \ e i(s(y)-s(x)) _ i| < 2. Therefore, (|2TT7|) is valid for 
s = 1. Furthermore, if we assume the validity of (|2.17p for s — k then since 
(2-18) 

a\ +1 {x) - a\{x) = % r V (y)^[e^- S ^ ~ l}V 2 (y)(a k ( y ) ~ a k -\y))dy, 



we have 

\a1 +1 (x) -a k (x)\ < ^ f ° \V (y)^V 2 (y)\p k (y)dy 
(2-19) l x 



h f X ° 

= ~~ k lJ p'(v)p k iy) d v 



P k+i (x)h k 



(k + iy. ' 

We would like to have an estimate for the function p(x) , therefore it suffices that we 
obtain an estimate for \ Vq (y) | ~ \ V 2 (y) \ dy. As already noted above, for x > x$ 



'x 



p(x)< \V Q (y)\^\V 2 (y)\dy<C(p) < oo 



■fi) 



Hence, the solution to the integral equation (|2.14p is given by 

oo 

(2.20) ai {x) = J2( a i( x ) ~ a t\ x )) 

s=l 

and satisfies sup x>Xo |ai(x)| < C(p) < oo uniformly in < h < 1. To derive the 
estimate for a[ (x) , we observe that 

l e US(y)-S(x)) = l V( ^(y) dye mv)-S(.*)). 

Therefore, using this observation and integrating by parts in (|2.10|) yields 
- \ £ {^py^ iSiv) - Six)) l]Vo(y)^V 2 (y)dy 

+ jTletW^W) - 1] (|jgj) V 2 (y)a l(2/ ) dy 
= -\ f V (x)i[ei^- S ^ ~ 1] dylVoiyy'Viiy)} dy 

_l [e i w ^(.))_ 1]( w)* %(afl) 

+ f X [e^- S ^ -l}(^) h V 2 ( y ) ai (y)dy. 



■V (y)< 

At this point we note that for x > xq (xq large enough), we have \d x Vo(x, ti)\ < 
Ci,^x ~ 2 ~ £ , uniformly in h. Hence 

(2-21) \dy[V {y)- l V 2 {y)]\<y-\ 

and using the boundedness of oi, together with that of p{x), obviously implies that 
for x > xq 

\a[(x)\ < C^aT 1 , 



ZERO ENERGY SCATTERING 



9 



uniformly in < fi < 1 as desired. For the case of a 2 one proceeds in essentially the 
same way using, however, the forward Green function rather than the backward 
one. This yields 

a 2 (x) = -J V Q (y)^[e^ s ^' s ^ - l]V 2 (y)[l + ha 2 (y)]dy. 

The same arguments as before now show that a 2 satisfies (|2.4[) . The Wronskian 
W{ip\ 1 tp 2 ) is obtain by evaluating at x — oo. □ 

The same analysis of course yields zero energy solutions with the correct asymp- 
totic behavior as x — > — oo. Note that the solution ip 2 is the sub-ordinate one, i.e., 
it is the unique (up to a nonzero scalar multiple) solution with the slowest possible 
growth. Hence, a zero energy resonance in the context of Theorem [1] would mean 
the existence of a nonzero solution to P(x,hD)f — with f(x) ~ c±|:r|2~~ Q ± as 
x — > ±oo and with 

«± = \J\+P± h ~ 2 

It is easy to see that such a solution cannot exist if V > 0. Indeed, let ^ be a 
standard cut-off function with x(0) = 1 and set xr( x ) '■— x( x /R)- If a globally 
subordinate solution f(x) did exist, then 



= hmsup(P(-, HD)f, X r!) = / [(/')>) + V(x)f 2 (x)} dx 

R — >oo J 

implies that / = 0, which is a contradiction. 

3. The Liouville-Green transform for small energies 
In this section, we consider the equation 
(3.1) - h 2 fl{x) + V(x)f ± (x) = Ef ± (x) 

where V is as in Theorem Q] As explained in the introduction, we will use the 
Liouville-Green transform to reduce (|3. If) to a perturbed Airy equation. We begin 
with a statement of the formal aspects (i.e., not involving estimates) of this trans- 
form, cf. Chapter 6 in [TH] and Langer's papers [H]-[I1]. Henceforth, V, Vo are as 
in Theorem [T] Throughout this section x > 0. 

Lemma 3. There exists Eq — Eq(V) > so that for all < E < Eq one has the 
following properties: the equation Vq(x; K) — E = has a unique (simple) solution 
on x > which we denote by x\ — x\{E; K). With Qq :=Vq — E 



(3.2) C = C(x, E- K) := sign(x - Xl (E; h)) - / y/\Q (u,E;h)\du 



defines a smooth change of variables x i— > C for all X > 0. Let q := — Op-. Then 

- = C = Vv> and 

-h 2 f" + (V-E)f = 



V > °> 3§ = C = v 7 ?' and 



transforms into 

(3.3) -h 2 w(() - (( + h 2 V((,E;hj)w(0 
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under w = \/C f — Here ' — ^ and 

Tr 1 -i i \-2 ~±d 2 qi 

Proof. Let E > be such that Vo(x; H) = E has a unique pair of solutions denoted 
by X2(E; H) < < ft). It is clear that (|3.2p defines a smooth map away from 

x = xi(E;h). Taylor-expanding Qo(x, E; K) in a neighborhood of that point and 
using that Vq(xi(E; h)) < implies that £(x, E; K) is smooth around x — x\ as well 
with C( x ii E;h)>0. Next, one checks that 

and thus, using -h 2 f" = {E - V)f, 



-h 2 w = q- 1 (E~V)w-h 2 q-- d2q 



d( 2 

q- 1 hQo + h 2 (x)- 2 /4)w-h 2 q-^w 

d 2 qi 



= ( w ( C )+h 2 (q- 1 (x)- 2 /4-q-^)w 

as claimed. □ 

We now analyze the properties of the change of variables introduced in the 
previous lemma. Recall that 

V (x; K) = (/4 + h 2 /A)x- 2 {\ + O(x^)) 

which implies that 



(3.4) x 1 (E;h) = c{h)E-i{l + 0{E^)), c{%) = yj p? + + h 2 /4 

It will be convenient for us to normalize the constants here so that c{K) = \J~2 
and we shall assume that for the remainder of this section. Moreover, we shall 
mostly suppress the harmless h dependence of various functions in our notation. 
We begin with the following "normal form" lemma which will allow us to describe 
the function £ in any region of the form x > eE~? (which, in particular, contains 
the turning point). Note that Lemma U normalizes the turning point xi to £ = 1 
by scaling out the energy. 

Lemma 4. Let e > but fixed. There exists a smooth map £ = E) on (y, E) € 
(e, oo) x (0, i?o) with Eq small so that £(E~z xi(E) , E) — 1 and for all (y,E) in this 
range, 

(3.5) l-E- 1 V {E-hy)=(^\l-C 2 ) 

and such that, with some constant £o(E), 

(3.6) ay,E) = y + Zo(E)+y- 1 p (y,E) 
(3-7) \d k E dlt(y,E)\<C k4 E- k y 1 -" 

for all k,i > 0. The functions £o ond po from (|3.6p satisfy 

(3.8) \dko(E)\ < C k E-\ \d k E dip (y,E)\ < C jk y^E- k 
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for all k,j>0 and uniformly in (y,E) G (l,oo) x (0,-Eo). For fixed < E < Eo 
the map y i— > £(y, E) is a global diffeomorphism whose inverse y — y(£, E) satisfies 
the bounds 

(3.9) y(^E)=Z + y Q (E)+C 1 MtE) 

(3.10) \d%dly(Z,E)\<C k4 E- k e- i 

for all k,£ > cmd with functions yo, po satisfying (|3.8[) 6u£ relative to £ rather 
than y. All constants are allowed to depend on e > 0. 

Proof. Set yi = y x {E) := E%xi(E). Then Vl = s/2 + 0(E?) as E 0+. Note 
that the O(-) term here satisfies 

d%0(Ei) = 0{E?- k ) Vfc>0 

due to the corresponding assumption on the error term of V. The same comment 
applies to every O(-) term appearing in this proof, both with respect to derivatives 
in E and spatial variable^. Define the change of variables £ = E) via 



(3.H) 
(3.12) 



f V sjl - E-WoiE-iu) du = t yjl-t- 2 dt y > yi 
Jyi Ji 

J V \J E^VoiE-iv.) -ldu = J Vt- 2 - 1 dt e < y < y x 



By monotonicity, these identities define a unique correspondence between £ and y 
on these ranges which is, moreover, smooth and strictly increasing on e < y < y\ 
and y\ < y < oo. By inspection, they also satisfy (|3.5|) . Since 

2 

it follows furthermore that the interval s < y < oo is transformed into one of the 
form < £i(E) < £ < oo where 



1 - E- x Vo(E-2 U ) = 1 - — + 0(Ei U - 3 ) 



(3.13) ^(E) = ^(O) + 0(E?) as £ -» 0+ 

and £i(0) > is a constant. We first show that the map £ = (;(y,E) so defined, is 
smooth for all (y, E) G (e, 2) x (0, E ) together with the desired estimates. To this 
end, write 

1 - E- l V (E-iy) = (y - Vl )U(y, E) 

U(y,E):=-E--z f 7 '(iH(yi + t(y - Vl ))) dt 



Then for all < E < E a and all k,£>0, 

(3.14) max \d k E d' y U(y,E)\ < C k , e E- k , min My,E) > c > 

l<a<2 " l<y<2 

For all e < y < 2 we rewrite p. lip and (|3.12[) in the form 
(3-15) (y- yi )Y(y,E) = (Z-l)X(Z) 



^We will say that a O(-) term behaves like a symbol if its derivatives are governed by such 
power-laws. 
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where 

Y(y, E):=U — t)U(yx + t(y — yi),E) dtj ' 



rl ^(2 + 5(^-1)) 

By the preceding, 



yl — v 2 dv 



\d k E d e y Y(y,E)\ < C Kl E-\ |^'X(0| < C, Vk,£,j>0 

uniformly on the interval e < y < 2 and the corresponding interval in £. By the 
inverse function theorem, (|3 . 1 5(1 defines a (unique) smooth map also locally around 
f = 1 and y = y\) this agrees with the previous definition for y ^ y\ and thus 
furnishes the desired smooth extension through the point y = y%. Furthermore, 
from 

yi = V2 + 0(Ei), d k EVl = 0(E^ k ) 
and 1(335]) . (|3~14|) we conclude that 

max 9 |5|^^( 2 /, J B)|<C w ^- fe 

l<y<2 " 

for all fc,^ > 0, < E < E . For large y, we write 

£{1-^(1+0^.-))}* 

The integral on the right-hand side satisfies 

K , 

VI -«- 2 dv = £ + K + oic 1 ) 

with a constant n, whereas the one on the left-hand side is equal to 

y + y a + 0{y- l ) + 0{E^-y- 2 ) 
with a constant y (E). It is easy to see that 

£ + k + OiC 1 ) =y + yo + Oiy- 1 ) + 0(E^y~ 2 ) 
implies (|3.6p and we are done. The statements about the inverse follow easily. □ 

We refer to Lemma 0] as a "normal form" since (|3.1|1 , on the interval x > E~~ i , 
turns into a suitably normalized (perturbed) Bessel equation in the variable £, see 
the appendix. By means of the change of variables introduced in Lemma [4] it is 
now an easy matter to describe £(x,E) from Lemma [3] on the interval x > E~z. 

Corollary 5. For all < E < Eq the following holds: there exists a constant 
c > \/2 so that on the interval e < \f~Ex < Co 

(3.16) ((x, E) = 2i (£-1) [l + 0(£-l)] 
with O(-) analytic and £ = £(y/~Ex,E). For all x > E~zcq, 

(3.17) l<;Hx,E) = z+ 7 +o(c 1 ) 

where 7 is some constant and O(-) is analytic. Neither of the O(-) terms here 
depend on E (other than through £). 
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Proof. We begin with £ close to £ = 1. The action S(x, E) then satisfies, with £ as 
in Lemma IU 

S(x,E) = aga(x - xx{E)) f ^/\E-Q {u)\du 

p\fEx i 

= sign(z - x x {E)) / V|l - E-W {E-iy)\ dy 

rf(y/Ex,E) 

= sign(Z(VEx) - 1) J V|1-^ 2 I^ 
= ^sign(£-l)|£-l|i(l + 0(£-l)) 

where the O(-) term is analytic in |£ — 1| < 1 and £ = ^(^/Ex, E). In terms of £ 
this means that 

C(as,£)=2*(£-l)[l + 0(£-l)] 

which is (|3.16|) . The constant cq is chosen so that 1 < £(cqE~2,E) < 2 for all 
< E < E . Since xi(E) = \/2/.E + o(l) as £ ^ we see that c > V%. As for 

HOT, 

£,(VEx,E) 



S{x,E) = I y/l - T 2 dv 

rt 

(l + 0(7r 2 ))rfr, = e + 7 + 0(r 1 ) 

Ji 

Since C = (§£)^ 5 we are done. □ 

In the region < x < exi(E) we have the following description of £(x, E) with e 
the same as in Lemma SJ In fact, in the following lemma we will need e small and 
then use this choice in Lemma |4j 

Lemma 6. For sufficiently small and fixed e > there exists a smooth function 
x(x, E) on < x < exi(E) with 

x{x,E) = x(l + 0(Ex 2 )) 

and such that 

2 r x ^i E ) 
(3.18) ~(*(x,E)= y/V (v) dv + 0{E log E) 

Jx(x,E) 

for all < x < exi(E) and < E < Eq. The O(-) here behave like symbols. 
Proof. Define x via 

rx{x,E) 



y/V (v) dv = I y / V (u)-Edu = J y/V a (u)(l + 0(Eu 2 ))du 
y/V (u) du + 0{Ex 2 ) 



Provided e > is sufficiently small (independently of E, of course), it follows from 
monotonicity considerations that x(x,E) exists with the desired properties. Next, 
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note that for all < E < E , 



I 



x x {E) px 1 (E 



y/V {u) du- / VWu) -Edu = 0(E log E) 



and thus 



2 



xi{E) 



-C 3 =/ y/V (u) -Edu- / y/V (u)-Edu 
Jo Jo 

xi(E) rx(x.E) 

y/V (u) du - / y/V {v) dv + 0(E log E) 
Jo 

xi(E) 



y/V {v) dv + 0{E log E) 

x(x,E) 

as claimed. □ 

The point of (|3. 18|) is that 0(E log E) is negligible as compared to the integral 
on the right-hand side which is on the order of | log(E(x} 2 )\ > 1. Thus, £ behaves 
to leading order like | log(E(x) 2 )\^ . We now turn to estimating the functions q, V 
from Lemma El In what follows, the notation A ~ B will denote proportionality of 
A, B > by some constants that are only allowed to depend on V. Also, A < B 
will denote A < CB where C is a constant, and similarly for A > B. 

Lemma 7. Using the notations of Lemma\3[ let < E < Eq. Then on the interval 
C > — 1 f/ie functions q = g(£, £7) and V = V(£, U) satisfy 

\d k E dlq\<C k4 E l - k {C)- 1 -" 

|d|^y(c,£)| <c MJ B- fc (c)- 2 -^ vm>o 

On i/ie interval £(0,E) < £ < —1 we wezw g, as functions of x via (|3.2[) . TTien 
one /ias g ~ |C| ( x )~ 2 an d there is the representation 

(3.20) V(C, £0 - -y^2 + ^{EPoix, E) + (z)- 3 /^, £)) 

where j3j satisfy the bounds 

\d k E diP 3 (x,E)\<C k , e E- k (x)- e j = 0,l 

\d k E diq\<c k/ E- k \c\- l {x)- 2 -" Vk,l>0 

All constants are independent of E. 

Proof. The case C > — 1 corresponds to x > exi(E) by Lemma 2] and Corollary [5] 
We now use that corollary to write 

C = Z = Z(y,E), y^E^x 

Then 

q = (C) 2 = E{dtat, E)) 2 (d y ay, E)f ~ £ 

The derivative bounds on q now follow from those obtained in Lemma [Hand Corol- 

jL 



laryO As for V, we compute, with ' — " 



(3-22) V = ia-^)- 2 - fl"*^ = \<l- l {*)-* + ^TV - \<rH 
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From the bounds on q which we just derived, the last two terms on the right-hand 
side of (I3.22| are < £ ~ 2 and behave as stated under differentiation. To treat the 
first term, we invoke (|3.9[) . (|3.f 6|) , and (|3.f 7p to write 

q-'ix)- 2 = q-^E-iy)- 2 = q-^E^y^, E))~ 2 

= q- 1 (E-hm,E),E))- 2 

which implies the correct bounds. Indeed, if |C| < 1, then q ~ E and the change 
of variables £ i— ► £ i— > y has derivatives of size < 1 relative to £ uniformly in E. 
This implies that q _1 (a;) _2 ~ 1 with derivatives with respect to £ of size < 1; 
furthermore, each derivative in E costs one power of E. Next, if £ > 1, then the 
change of variables ( i— ► £ i— > y acts like C 3 by Corollary [5] Thus, 

with each £ derivative gaining one more power of decay in (. 

In the remaining case ^ < — 1 one first calculates, on the one hand, 



. d 2 q~< 



16C 2 
5 



l Qo _ l_Qo_ _ 6 /^oy 

4Q 8CQo 16Vq / 



where 



d_ 

dx 



16C 2 

Thus, from (|3~2"2")l . 
1 



9 



Qo 8(Qo 16 VQ, 
4Q 16VQ 



_i d 2 q~< 



dc, 2 



where we have selQ 
(3.23) 0o(x,E) 



16C 2 
5 

16C 2 " 



4Q 16VQo 



Ep {x,E) + (xy 3 f3 1 (x,E) 



(3.24) 



/3i(a:,f0 := ( x ) -;(*) 



Ulvb Q o y 16VVQ 7 Vv 

-4 W 4Fo 16 V Vb / J 

As already noted in Section [21 the x~ 2 terms inside the brackets in (13.24j) cancel 
so that the leading order is x~ 3 . In fact, \d l x f3i(x, E)\ < Ce(x)~ e in view of our 
assumptions on V, see Theorem [TJ As for (3q, we note that in the range ( < — 1, 
one has Qq ~ Vq. Since Qo = Vq — E, this implies that the expression in brackets 
in (|3.23|) is < E together with the natural derivative bounds. The bounds on 

q = (V - E)^- 1 ~ (x)- 2 \C\-\ |CI ~ I log(E(x) 2 )\i 
follow from (|3. 18|) and we are done. □ 

In Lemma[7]the modification of V to Vb only played a role in the regime £ <C — 1 
which is the same as x < ex\(E). This is natural, since we know from Section [2] that 
this modification really comes from the E = case which corresponds to x\ = +oo. 
We will see this mechanism at work in the following section, too. 



j9l does not depend on E, but this makes no difference. 
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4. Solving the perturbed Airy equation 



This section is devoted to solving (|3.3|) . at least in the asymptotic sense relative 
to h. We shall use the notations and results of the previous section. For the 
properties of the Airy functions Ai, Bi listed below we refer the reader to Chapter 11 

of mg. 



Proposition 8. Let ho > be small. A fundamental system of solutions to ()3.3j) 
in the range C < is given by 

cf >1 (<:,E,h)=Ai(T)[l + ha 1 (C,E,h)] 

MC,E,h)=Bi(T)[l + ha 2 (C,E,h)} 

with t := —h~i(,. Here a^o^ are smooth, real-valued, and they satisfy the bounds, 
for all k > and j = 1, 2, and with Co := £(0, E), 

\d% aj ^, E,h)\< E' k min [ft* (ft" *£) a , l] 

\d k E d caj {Q,E,K)\ <^- fc [R-»<fi-'C)-*X[-i<C<o] + KI*XKo<C<-i]" 
uniformly in the parameters < h < ho, < E < Eq. 

Proof. Let 01,0 (C: ft) := Ai(r) an d 02,o(C, ft) := Bi(r). We seek a basis of the form 

fc(0 = <KC, ft, £) = &,o(C, + ftoi(C, ft, £)) 

for £ < This representation is meaningless for C > since 0j,o have real zeros 
there. On the other hand, on C < they do not vanish. We obtain the equation 

1: 



(4.1) 



(4.2) 



Oj,0«j) 



-V<t>io(l + haj) 



for j = 1,2 where ' = 8q. A solution of (|4. 2|) on C < is given by, with 02(C) 
02(C) ft, #)> 



02(C) == - 



(4.3) = -m 



<t>\fi{n,ti) / &2l{Vift)d-qV(T 1 ,E){l + ha2{ri))dTi 



-h~3(, 



Bi 2 (it) 



Bi _2 (t;)dt; V(-fi*u, £)(! + fto 2 (fi 5 u)) du 



This solution is unique with the property that 02(0) = o 2 (0) = 0. Recall the 
asymptotic behavior, see [Hi] . 



Bi(x) = 7r *x l e z x2 [l + 0(x 2)] as 1 
Bi(x) > Bi(0) > V x > 

1 1 1 2 3 3 

Ai(ar) = -7r~ saT^e - ^' [l + 0{xT^ )] as a; 
Ai(ar) > V x > 
Also note the useful fact, valid for any < xo < x\, 

fV 2 ( y )d y = W^-^4 

J XQ \Bi(x ) Bi(xi) 



00 



00 



(4.4) 

which implies that 



Bi 2 (x ) f 1 Bi' 2 (y)dy 

J in 



< 
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The leading term in (|4.3|) . i.e., 

i /— 

a 2 , {C,E,h) := -h* / Bi 2 (u) 



Bi~ 2 (v)dv V{-hzu,E)du 



therefore satisfies the bound (dropping E, H from 02.0 for simplicity) 

2 



|02,0(C)| <^ 



(u)~3|y(_^3 U;J B)|dw 



We now use the estimates from Lemma [7] to bound the right-hand side. If — 1 < 
£ < 0, then this yields 



(4.5) 



MO! 



On the other hand, if £0 : = C(0) E) < ( < —1, then we obtain 



(4.6) |a 2j0 (C)l < ti3 



(4.7) 



(u) 2 du 

2 



a 



n 3 



du 



The variable z appearing in (|4. 7[) is tied to the integration variable u via —hiu~ 
C(z,E), see Lemma [3] The integral in (I4.6[) and the first term inside the brackets 
in (|4.7[) contribute 



/l~3 



-n.-3( 



a ' 2 {h^u)~ 2 du < 1 



?i~3 



Next, with C(^2, £■) = — 1) and v = C,{z 1 E), dv = ^Jqdz 



h 3 



-C 



u 2 q L (h~3u)Edu = E I v 2 q 1 (v)dv 



E 



>(z,E)dz <E zdz 



<Exi(E) 2 < 1 



Finally, using that 4^ = \/Q once again one obtains 



Hz I u 2 q 1 (h^u)z 3 du — v 2 q 1 {v)(z) 3 dv 



< 



Q *(z,E)(zr 3 dz 
(z)- 2 dz < 1 



In summarjQ, 

\a 2 AC,E,h)\ <min(l,7P</r§C>') 

^Had we used V instead of Vo in our definition of then we would be losing a factor of logi? 
at this point. Indeed, for the case of V we would need to replace E+ (z)~ 3 by the strictly weaker 
{z)~ 2 in 114.71 1 which then leads to the logarithmically divergent integral f^ 1 (z)~ 1 dz. 
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uniformly in ( 6 [Co , 0] , < E < Eq , and < H < Hq. Due to the linear nature 
of (|4.3|) . a contraction argument now yields the same bound for a 2 ', in fact, due to 
the derivative bounds of Lemma we obtain the more general estimate 

\d%a 2 ((,E,h)\ <C k E- k mm(l,hi(h-i()^) Vfc>0 

uniformly in the parameters. As for the first derivative in £, observe that 

(4.8) da(C) - -T^TT T <t>lo(v, E)(l + ha 2 ( V )) d V 



whence, for all — 1 < £ < 0, 

2 



<h-^Br 2 {-h-^c) f Bi 2 

Jo 



(u) du 







ft -1 ICI* 



<ft-3(?l-IC)5e-3 fi ^l* I "' (v)-$\v\-iei v dv 

Jo 

<h-^h-i&e-*' l - 1 M i {h- 1 \<;\i)-h*' l - 1 M i 



If Co < C < -1) then 

|d 2 (C,£,7i)|</r^o(C,ft) I ^, (v,K)\v(v,E)\d v 



c 



(4.9) < h-*Br 2 (-nr%() / Bi 2 (u)du 

Jo 



h 3 C _ r „ E- 



r~ * r 
(4.10) +fi-3Br 2 (-ft-lc) / 2 Bi 2 (w) (0u 



du 



where z has the same meaning as in (|4.7[) . First, note that (|4. 9[) is rapidly (in fact, 

3 

super exponentially) decreasing as |£| increases: I gill <e~KI 2 . Second, we bound 
the first part of (j4~T0l) by 

r-ft 3 c 

fr*Br 2 (-?rzQ / 2 Bi 2 (u)(fr§u)- 2 du 

2 



< h-?Br 2 (-h-i() [ u-h? ul du 
<?H(/rl|C|)-*<|C|-' 
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The contribution to (14.101) involving E is 

h-?Br 2 (-h-?() / 2 Bi 2 (^(qihiu^Edu 

_ 2 

3 r-h 3C 3 

(4.11) <Eh-^{h-i\C,\)^e-^^ / 2 ti-'e^ (g(^u)) _1 du 

(4.12) < J E/r 1 |C|*e-& lcl * ^ C ze^ t ' i dz < |C|^(a;) 2 

To pass from (|4. 1 1|) to (|4.12|) . we substituted u = h~%v and then changed variables 
cfo = \J dz followed by vq(—v) = Qo(z); in (|4.12[) the relation between v 
and z, as well as £ and x, is given by (|3.2[) . i..e, u = £(z, £ = £(x, To pass 
to the final inequality in (|4.12[) we integrate by parts so as to gain a factor of ft: 

J~ C z e^ vi dz < £ z e« £ 1 VO^) *» dz 

< h{x) 2 e^^ y/QM)*t = h(x) 2 e^ 

where £(^2, ft) = —1. Finally, we turn the contribution of (z)~ 3 in (|4.10p . Using 
the same conventions regarding the relation between the variables this contribution 
is of the form 

_ 2 

ft" 3 BP 2 (-ft" I C) / 2 Bi 2 (M)(g(ftSw))" 1 (z)- 3 du 
<ft-3(ft-i|C|)3e-AlCI § / 2 it-sel" 1 ^^!^))- 1 ^)^^ 

< ft^lCl^e-AKI 1 ^ f ( z )- 2 e 37T^ rf z < |C|^(x)- 2 
The final inequality here is based on the same kind of integration by parts as before: 
r\z)- 2 dz < [ X \zT 2 el f** 1 VQ^M *j dz 

Jl Jx 

< h{x)- 2 ei S* 1 VQ^) d » = h(x)- 2 e^ 

with X2 as above. In conclusion, we estimate the contributions of (14. 9|) and (|4. 10|) 

by 

(4.13) \a 2 ((,E,h)\ <fi-'(ft-'C)"'X[-i<c<o] + lCl'X[Co<C<-i] 
as claimed. 

Next, we turn to (j>i((^,E) (dropping ft for simplicity). As usual we make the 
reduction ansatz 

MC,E)=g(t,E)Mt,E) 

which leads to the equation (<f> 2 g) — 0. At this point it is convenient to extend the 
solutions <j>2, which are originally defined on the interval £(0, E) < £ < 0, to all of 
£ < 0. This is done in such a way that the bounds (|4.1[) remain valid for £ < £ 
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without, however, making any reference to the ODE (|3.3j) for those £. We can now 
solve for for g in the form 

f c 

J — oo 

Inserting our representation of fa into this formula yields 

rC 

fa{(,E)=Trh-%Bi(-h-%()[l + ha 2 {(,E)} / Bi~ 2 (-H^ r,)[l + tm 2 {r,,E)]- 2 dr] 

J —oo 



First, we note that from 



7Tfi-3Bi(-7r3£) / Br 2 {-tr^7i)dri = Ai(-hr*Q 

J —oo 

Second, [1 + ha 2 ]~ 2 = 1 + Ka 2 where a 2 satisfies the same bounds as a 2 (since 
| d2 1 ^ !)■ Thus, inspection of our formula for fa reveals that a,\ = Tr(a 2 + ai) where 

oi(0:=fi -S |j(-ft-»C)[l+fr»2(C.-E)] f Br 2 (~h-ir 1 )a 2 (v,E)di 1 

_ Bi 
~ Ai' 



{-h-H)[l + ha 2 ((,E)} / Bi-^^aa^??,^)^ 



ft~3£ 



Furthermore, from (|4.4j) . 

9 Br 2 (r/)S 2 (^7 7 ,^) ( i7 ? = - 



'-ft~3f 

(4.14) =^(-^iC)S 2 (C, J B)-ft | 



2 r Ai 

, a 2 {-h*<n,E)d —(rj) 

^(r?)(a 1 a 2 )(-?l3?7,£;) 
Bi 



where <9i refers to the derivative in the first variable. The first term in (|4. 14[) makes 
an admissible contribution to a\ whereas the second one is controlled as follows: 



Ai 

*-§C Bi 



(rj)\(dia 2 ) (-hsj],E) I dr] 



<hh^ 



K~3( 



< h^{h^^()^ 1 X[-i<c<a] +%[Co<C<-i] < ^(^C)^X[-i<C<o] +X[Co<C<-i] 
as desired. For the derivative in £, 



/oo 
2 Bi 
-ft - 3 C 



i 2 (r))a 2 (h3i],E)dr] 



+ h-%(AiBi)- 1 (-h-ic)[± + fia2((,E)}a 2 {(,E) 
+ h J-(- h ~H)a2(C,E) 



Bi 2 (r))a 2 (h^r],E)dr] 
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Using (|4. 14|) we remove the dangerous ftri terms from the first two lines here 
whence 

/OO A ■ 

-ft - ' sC ai 

Bi f 00 
(4.15) +n— (-n-tC)ia(C,^0 / , Bi^^ffi^.E)^ 

The contribution by the first integral here is treated as the integral in (|4 . 14[) and 
is bounded by 

~ fr^{%-i C>~^X[-i<c<o] + ICI*X[Co<c<-i] 

which is exactly as needed. Finally, the contribution of (14. 15[) is bounded by 



< h 



X[-i<C<o] + lCl 2 X[Co<C<-i] 



~H0 2 ^ h 3 ( h 3 C> 2 X[-i<c<o] + ICI 2 X[Co<c<-i] 

and we are done with the k = case of (|4.ip for a\. However, since E enters into 
a\ only through 02,02 which do satisfy (|4.ip for all k > 0, we see that the previous 
estimates carry over unchanged and provide the stated estimates for d E ai((, E) 
and d%d cai ((,E). ' □ 

We remark that the method employed in the previous proof does not extend easily 
to derivatives d^aj with £ > 2 (that is, without losing excessive powers of ftr 1 ). In 
principle, it is possible to treat £ = 2 by a similar method, but instead of a sharp 
cut-off at £ = one needs to use a smooth cut-off function in (|4.3p . However, 
the calculations are quite involved and it is not clear how to extend this approach 
systematically to higher £ (the same comment applies to Proposition [9] below). On 
the other hand, for the purposes of Theorem [TJ as well as for those of [23] and [24], 
it suffices to treat the first derivative in £ (however, we do need many derivatives 
relative to E). Next, we turn to ( > which requires an oscillatory basis. 

Proposition 9. Let Hq > be small. In the range ( > a basis of solutions 
to (|3.3p is given by 

MC, E; h) = (Ai(r) + iBi(r))[l + fi6i(C, E; h)] 
MC, E; K) = (Ai(r) - iBi(r))[l + hb 2 ((, E; K)] 

with t := —h~i( and where 61,62 are smooth, complex-valued, and satisfy the 
bounds for all k >0, and j = 1,2 

, A1a , \d h E b 3 (CE;h)\<C k E- k (C)-l 

(4.16) , , , 

\d c d k E bj(c,E)\ < c k E- k h-*(h-H)-HQ- 2 

uniformly in the parameters < h < Hq, < E < Eq, C^O. 

Proof. Let Vi,o(C;^) := (Ai + iBi)(r) and ip 2 ,o(C,^) ■= (Ai - iBi)(r). We seek a 
basis of the form (dropping h as an independent variable from the notation) 

Vj(C) = E) = ^,o(C)(l + W>i(C, E)) 
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for £ > 0. This representation is meaningful since Ai and Bi have no common zeros 
(as their Wronskian does not vanish). We obtain the equation 

(4.17) (^>i) = -^o(l + fiij) 

for j = 1,2, c.f. (|4~2|) . A solution of ([4~T7| on C > is given by, with bj(() = 



(4.18) 



-1 



b 3 (0 ■■= -T %M / Wri, E)(l + hbM) dv 



Recall the asymptotic behavior, see [19] . 

(4.19) Ai(-a;) ± iBi(-x) = -JL^^ - *^! + 0(f -1 )) 



7T2 x 4 



as a; — > oo. Here £ = |a;2 and the O(-) term is complex-valued and exhibits symbol 
behavior: 



deO(C L ) = 0{C 1 ~ k ) Vfc>0 



Therefore, for any > x$ > x%, 



(Ai + zBi) 2 (a;i) / (Ai + zBiT 2 (y) dy < (xi)' 



The leading term in (|4.18|) , i.e., 



b lfi (C,h,E) := TP 



(Ai+iBi) 2 (-w) 



(Ai+iBi)~ 2 (-v)dv V(-hzu,E)du 



therefore satisfies the bound, see Lemma [71 



IMOI ^ h ~ 3 



< TP 



OO 
2 

fi"3 C 



(u) _ *|V(-7Pu,.E)|du 



uniformly in £ > 0, < E < Eq, and < h < ho. Due to the linear nature 
of (|4.18p . a contraction argument now yields the same bound for b±; in fact, due 
to the derivative bounds of Lemma [Jj relative to E, we obtain the more general 
estimate 

|d|MC,£)| <C k E- k (cr~ 2 Vfc>0 
uniformly in the parameters for both j = 1,2. As for the first derivative in £, 
observe that 

h- 1 



(4.20) 



&i(C) 



€o(0 



^ t0 (vMv, E^l + hb^dr] 



In order to exploit the cancellation in this integral, one integrates by parts once. 
To this end, write for u > 0, 

(Ai + zBi) 2 (-w) = e%M* lj(u), \w(u)\ < («)-', |w'(«)| < 



Since 



^,o(C;^ = ef<^>M?i-!c) 



^, (C; ^ rfC = ^ (ft-»C>" s w(ft-»0 d 

2i 
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integration by parts yields 

(4.21) =o(^(/Hc)-%r 2 )- 

(4.22) ~2^(C)/ C ^"^M^r^) -2 ^)^ 

The leading order here is given by (|4.2ip : indeed, if we estimate the 61(77) term 
in (|4~22|) by (l4~2Tj) . then (j4~22)) < ft(C)~ 4 , which is much better than (|4%|) . The 
conclusion is that 

|a c fe(C,^)l < s- fe ft-3(fi-§o-'(C>- 2 

as claimed. □ 

5. The proof of Theorem Q] 

Let f±(x, E; h) be the Jost solutions of P(x, HD) from (jl.ip . For ease of notation, 
we shall first assume the symmetry V(x) = V(—x) and later indicate how to treat 
the general case. Also, as usual, we drop h from the arguments of functions. Then 
f-(x, E) — /+(— x, E) so that the Wronskian of /+, /_ is 

W(E) = -2f + (0,E)f' + (0,E) 

Next, from (|3~19|) . and with ( = ((x, E) as in and T + (E) as in (fl~5| . 



This is obtained by matching the asymptotic behavior of /+ with that of ip2(() as 
x — > 00 and we used the relation w — f from Lemma [3 We now connect -02 to 
the basis <pj((, E) of Proposition [51 

MC, E) = Cl (E)MC, E) + c 2 (E)MC, E) 

where 

(n _ W(M;E)M,E)) (F] _ W(M;E),M;E)) 
11 ' W(M-:E),M-,E))' 2[ ' W(M;E),M;E)) 

By Proposition [5J 

w(M-iE),M-,E)) = -tr%W(Ai,Bi) + o(jrf) = -7r _1 7rt(i + 0(h)) 

where we evaluated the Wronskian on the left-hand side at ( — 0. Next, by Propo- 
sitions O and [51 

W(th('>E)> M, E )) = -h-i[(Ai(0) - iBi(0))Bi'(0) 

- (Ai'(O) - iBi'(0))Bi(0) + 0(h)] 
= -h-l[W(Ai,Bi) + 0(h)] 

W(M;E), <f>i(-,E)) = -hrl [(Ai(0) - zBi(0))Ai'(0) 

- (Ai'(O) - iBi'(0))Ai(0) + 0(h)] 
= -ft"! [iW(M, Bi) + 0(h)] 



T,(E) 
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so that 

(5.2) c 1 (E) = l + 0(h), c 2 (E) = -i + 0(%) 

where the O(-) terms satisfy \d^O(h)\ < Ck E~ k . For the remainder of the proof, 
we set Co := C(0,E). Thcn 

/+(0, E) = v^e^^+T) E*h~«q~* (Co ) "02 (Co , E) 



,,T + (E) . x x J 

V5Fe t (-^ + ^^ft-»g-»Kb)[ci(^i(Co,-B)+ca(£)0 2 (Co,^)] 



/;(0,s) = ^ z ^ £1 +5)^ft-k'(0)g-^Co)[V' 2 (Co,s)-^(Co)V'2(Co,s)] 



'ire 



l (— +4) £*ft-*C'(0)g-*(Co) [<*(£)#(&, + c 2 (^)^(Co, 



- i £ (C„) (d (£)& (Co , E) + c 2 (E)<j)2 (Co, E))] 

Recall from Lemma [3] that C' = q^ ■ From V'(0) = we obtain 

.,, v Qo(0) g(Co) 
Co Co 

and thus 

/+(0, £0/+(O. = ^ 5 e 21 — fi-3 [ Cl (^)^(Co, £) + c 2 (£)0 2 (Co, ^)]x 

x [ci(£;)0' 1 (Co,S) + c 2 (i?)^(Co,^)] 

+ | 7 r^e 2i ^ £1 S-iCo" 1 [ci(£J)^i(Co,£;) + c 2 (£)<MCo, E)] 2 

From Proposition [51 

^i(Co,^)=Ai(-a-iCo)(l + 0(R)) 
2 (Co,^)=Bi(-fi-iCo)(l + O(?i)) 

#(Co, £) = -S-tAi'(-r*Co)(l + 0(a)) + 0(»)|Co| 4 Ai(-fi-f Co) 

2 (Co,^) - -fi-»Bi / (-»-»Co)(i + 0(a)) + o(a)|Co|'Bi(-a-lco) 

which implies via the standard asymptotics of the Airy functions that 

1 (Co,S) = (4 7 r)-5( ?l -l| Co |)-l e -i^ 1 lCo|i (1 + o(n)) 

M(°,E) = h-H4n)-Hfr-hCo\) L * 6-^^(1 + 0(h)) 

<j> / 2 (C 0l E) = -h-^-Hh-^\C \)iei h ^a + O(h)) 

Hence, using that e'^ 1 ^ = 0(h) where d k E 0(h) = 0(E- k h), one concludes that 

h- 3 [ Cl (£)</»! (Co , £) + c 2 (£)0 2 (Co , E)] x 

x [ Cl (£)</,' 1 (Co,-B)+c 2 ( J B)^(Co )J B)] 

= 7r- 1 a- 1 ei ft_1 lf°l 5 (l + 0(a)) 
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as well as 

ft-ko" 1 [ci(i?)^i(Co,S) + C2(i?)02(Co,i?)] 2 = -7r- 1 |Cor i et R " 1 l^l i (l + O(n)) 
Since T(E) = 2T + (E) we finally arrive at 

W(E) = -2f+(0,E)f' + (Q,E) = -2ie 2lI± ^ E^h^e^ 1 ^ (1 + O(h)) 

(5.3) = -^le^^+^^il + 0(h)) 

We used here that 

^ICol 1 - 2 p ^V Q (rj)-Edr, = S{E) 

All the 0(h) appearing above behave as required under differentiation with respect 
to E; indeed, this is both due to the bounds of Propositions [5] and as well as the 
aforementioned fact that 



O(h\( \-?)=O(h) 



has the required behavior since 

\( \-2 = 0(\logE\-*) asE^0+. In view of (fl~ 
(|5.3[) implies the sought after asymptotic relation for Sn in Theorem [TJ see (jl.6|l . 
In order to find §12, and §21 (i-e., the reflection coefficients), we need to also 
asymptotically evaluate the following Wronskians: 



W(f+(; E), /_(•, E)) = W(f + (;E), f-(;E)) = -2 Re [/+(0, E)f + (Q, E)]. 
Using the same notations as in the computation of W(E), we obtain 
-2Re [/+(0, E)f' + (Q,E)] = -2Re {t^/H [ Cl (E)MCo, E) + c 2 (E)fo(Co, E)]x 



x[c 1 (E)cj ) ' 1 (( a ,E) + c 2 (E)<f ) ' 2 ((o,E)]} 
Finally, evaluating this expression as above, we obtain 



^Eih-Ho 1 \ci(E)MCo,E) + c 2 (E)MCo,E)\ 2 



2y/~E siMl 



W(f + (;E), f-(;E)) = -2Re [/ + (0, E)f' + (0, E)] = ~j~e^(l + 0(h)). 

Forming the ratio between this formula and the one for W(E) yields the desired 
expression for §12 = §21, see (|1.6[) . Indeed, 



r (E) = _ W(M;E),f-(;E)) 0(H)) 
[ ' W(f+(;E)J_(.,E)) U + ^jj 

where 0(h) behaves like a symbol with respect to E, as usual. This concludes the 
proof of Theorem [1] in the symmetric case. If V(x) ^ V(—x), then only minor 
changes are needed. Indeed, on x < we can still use the same bases <t>j,ipj 
from Section [4] but with £ = £(— x,E). This is due to the fact that the difference 
between the left-hand and right-hand branches of V does not affect the estimates 
from Section [4] (since we are assuming inverse square decay at both ends and the 
constants [i± have no effect on the leading order behavior). Let 

Co(-B) f :=§/ VVo(v)-Ed V 

L Jx 2 (E) 
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Thus, in addition to the expressions for f + (0,E) and f' + (0,E) from above we now 
also have 

/_(0,25) = \/^e i( T+f) Eih-iq-i(Co)[ci{E)(f>x(Co,E) + C2(E)MCo,E)} 
/i(0,25) = -V^^+S) Ein-k'(0)q-^(Co)[ci(E)<P[(C o ,E) + c 2 (E)cp' 2 (Co,E) 
-^((o)(c 1 (E)^ ,E)+c 2 (E)M(o,E))] 

Inserting these expressions into 

W(E) = f + (0, E)f'_ (0, E) - f' + (0, 25) /_ (0, E) , 

and using that 

o [Co 2 + Co 2 ] = / VVo(v)-Ed v = 5(25) 

as well as T(E) — T + (E) + T_(i5), one again arrives at (|5 ,3[) . The same comments 
apply to the off-diagonal terms of the scattering matrix and we are done. As for 
the very last claim of the theorem concerning Vb = V + h 2 Vi, simply note that the 
main calculations entering into the above proof only make use of the leading order 
part of Vi, i.e., \{x)~ 2 whereas the cubic piece gets absorbed into the error term. 

6. From small to large energies 

In this section, we present an extension of Theorem[T]to the case of large energies. 
More specifically, suppose V is as in Theorem Q] but with the following additional 
properties: 

• < V(x) < 1 for all iel, V(0) = 1, V'(0) = 0, V"{0) = -1 

• V'{x) < for all x > 0, V'(x) > for all x < 

Note that this is precisely the kind of barrier potential considered by Ramond [22] 
(but without any analyticity assumptions). For the purposes of this section we refer 
to it as a simple barrier potential. Even though Theorem[T]by design only considers 
small energies < E < Eq, it is natural to ask to what extent it remains correct as 
Eq — > 1. As already remarked before, for energies 25 > e > there is no difference 
between V and Vb as far as Theorem [T] is concerned. Indeed, switching from V to 
Vb only affects the error term. Moreover, for the kind of V we are considering here, 
the theorem remains valid in any range < E < 1 — e with e fixed. This is due 
to the fact that in this range there is a unique pair of turning points x 2 {E), x\(E) 
as before. The action 5(25; K) lies between two positive constants (depending on e) 
and the previous proof goes through without changes. Somewhat more interesting 
and very relevant for later applications, cf. [33], [21], is the case where e = h a . The 
question is then how large a > can be allowed to be. First note that we can no 
longer expect the error term in (|1.6[) to be of the form 0(H) in that case. Rather, 
it will need to be 0(h s ) for some S — 5(a) > and this condition will determine 
how large we can take a. It turns out that the range < a < 1 is admissible here. 
In the following corollary, we use the notations introduced in Theorem [T] 
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Corollary 10. Let V be a simple barrier potential. For every < a < 1 there 
exists and Hq = Ho(a) small such that for all < h < ho and < E < 1 — K a 

Su(E; h) = e -mm)+iT(E;H))^ + h{l _ E) -l aii{E . 

' ' S 12 (E; K) = -te-^ T +^ h \l + - E)- 1 a 12 {E; ft)) 

and the correction terms satisfy the bounds 

(6.2) \d k E a ll {E-h)\ + \d k E a l2 {E-K)\<C k ^{E'\{l-E)- k ) V k > 0, 
with a constant Ck that only depends on k and V . 

Proof. We will only sketch the proof as there is no point in repeating all the details 
of the proof of Theorem[TJ In fact, inspection of the previous section shows that the 
main issue is to prove that Propositions [8J and [9] remain valid albeit with errors of 
the form h}~ a rather than h (we need to pay particular attention to the derivative 
8q) when E = 1 — h a . We will freely use the notation from Section [3] and 2J By 
the preceding comments, it will suffice to consider the range 1 — e<E<l — h a . 
In fact, it will be enough to set E = 1 — H a so that x\(E) ~ ft"?. The range 
< x < x\(E) then corresponds to the region — ft~3" < Q < 0. A simple calculation 
shows that q ~ h% in that range, as well as |V| < hr^ with the usual behavior 
under differentiation in E. In fact, for all < x < x\(E) we have 

V(x) -E= - I V'(y) dy ~ x\ - x 2 ~ (xi - x)xi 

J X 

and thus 

/• h ^ V-E xi(xi-x) § c, 

C~-xf(xi-x), q= — ~ — =xf~ft3 

' x| (xi — x) 

as claimed. Next, recall (I3.22p . viz. 

(6-3) V = \q~Hx)-* + A 5 -V - \q~H 

Since q = q~zq' where q' — ~ Xj 3 ~ , the second term here is of size 

TV < <T 3 (<z') 2 < T 4 ~ S"* 

The other two terms are smaller whence |V| < ft - ^ as claimed. Turning to Propo- 
sition [51 we seek a basis of the form 

ct> 1 (C,E,h)=Ai(r)[l + h 5 a 1 (C,E,h)] 

MC,E,h) =Bi(r)[l + h s a 2 ((,E,h)} 

with <5 := 1 — a. Proceeding as in the proof of Theorem[TJ we arrive at the following 
analogue of (|4.5[) 

i—nii 

|O2,0(C)l £ ^ / (u)-5|y(-^3 U , £7)1 rf u 



which yields 

Ko(C)l < hl- s (h-i()ih-^ < h 1 —- 5 

This shows that with our choice of 5, we have 

sup |a 2(0 (C)l ^ 1 

C(0,B)<C<0 
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For the derivatives, the analogue of (|4.8|) . viz., 

^ = / <Plo(v, K)V(ri, £0(1 + fotefo)) df 7 

yields 

2 

\a2((,E,h)\ < h*- 5 Br 2 (-h-i() / Bi 2 (u)7r^<frt 

./o 

<^|-*-¥{/i-lc)-5 <?r§ 

where we again used that a < 1 in the final step. An analogous estimate holds 
for 4>i . We claim that these bounds are sufficient provided the same type of estimates 
hold for the analogue of Proposition [9] at £ — 0. Indeed, inspection of (|5.1|) shows 
that in that case 

W(M;E),M-,E)) = -fi-i[(Ai(0)-iBi(0))Bi'(0) 

- (Ai'(0) - iBi'(0))Bi(0) + 0(h 6 )] 
= -ft-t[W(Ai,Bi) + 0(fi 5 )] 

W(M; E), fa (-,E)) = -H-i [(Ai(0) - ffii(0))Ai'(0) 

- (Ai'(0) - iBi'(0))Ai(0) + 0(h s )} 
= -»-t[iW(Ai,Bi) + 0(h 5 )} 

Note that there is an exact balance here between the h~i coming from the deriva- 
tives of the main contributions and the losses stemming from dj, bj. Hence, 

d(E) = 1 + 0(h 5 ), c 2 {E) =-i + o(h s ) 

as desired. Since fi. _1 |Co|^ ~ ft" -1 and thus also 

e ->r 1 |Ca|* _ Ofi 1 - 01 ) = 0(h s ) 

the reader will easily check that the remainder of the proof in Section[5]goes through. 

It therefore remains to deal with the oscillatory regime. In analogy with Propo- 
sition [9] we seek a basis 

MC.E; K) = (Ai(r)+iBi(r))[l +h%({,E;h)] 
MC, E- h) = (Ai(r) - *Bi(r))[l + H%((, E- h)] 

For this we need to understand V on £ > 0. First one checks that for all x > xi(E), 

{xf (x — x\) x\ < x < 2xi 
x% 2xi < x < 1 

13 X > 1 

and thus 

f £3 X\ < x <C 1 

L x 3 x > 1 

Hence, (|3.22p implies that 

\V\ ^ X ~^X[xy<x<\] + C 2 X[x>i] 
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Going through the proof of Proposition [9] shows that 

l^(o)| < h 1 -"- 5 < i, 1^(0)1 < 0- s -^ < JT* 

as desired. The derivatives relative to E are left to the reader. □ 

Thus, the semi-classical approximation obtained in Theorem [TJ breaks down pre- 
cisely at E = 1 — h. As is well-known, the Airy equation is no longer the correct 
approximating equation for energies close to the unique maximum 1^(0) — 1 of a 
simple barrier potential. In fact, there exists an analytic change of variables which 
reduces the Schrodinger equation with such energies to the Weber equation locally 
around the origin. Alternatively, Ramond [22| invokes micro-local methods and the 
Hclffcr-Sjostrand normal form in that case. 

Appendix A. A normal form reduction to Bessel's equation 

In this section we sketch an alternative route for the asymptotic analysis of 
Section [4] It is based on Lemma [4] and reduces equation (|3.1|) to a Bessel equation 
rather than an Airy equation. However, we emphasize that these approaches are in 
fact quite related as the Airy functions are used to describe Bessel functions J n and 
Y n in the large n — hr 1 asymptotics very much in the spirit of Section 21 see [19] . 
A possible advantage of working with the Bessel representation lies with the fact 
that they apply to all x S [eE~ 3 , oo) which is a region containing the turning point 
x\{E). On the other hand, since they cannot be used on the region [0,eE~z], one 
is again faced with a connection problem as in Section |4j Moreover, we have found 
that using distinct changes of variables in these two regions leads to a number of 
complications as compared to the global action-based coordinates introduced in 
Lemma El For this, as well as other reasons, we ultimately found it technically 
advantageous to work with the Airy approximation directly, but we wish to sketch 
the Bessel method since it seems to be of independent interest. In this section, we 
shall use the notations of Lemma [4] and always work on y > 1 which transforms 
into £ > £,i{E), see (|3.13p . First, a preliminary technical lemma. 

Lemma 11. The function /j,(rj,E) := (d^y^q, E)) 2 (d yy £,)(y(r], E), E) satisfies 

(A.l) \d k E d^{ Vl E)\<C k] E- k n-^ 

and the positive smooth function 

(A.2) Cl(r], E) := exp ( - J n{t,E)dt} 

satisfies ^- = fi and 51 = 1 + 0(rj~ 2 ) (we write ' = d n ) with a symbol-type 0(r]~ 2 ). 

Proof. The rj~ 3 decay in (|A.1|) is due to (|3 . 9[) . Otherwise, the lemma is an imme- 
diate consequence of Lemma |4j □ 

Now for the transformation of the equation with V, Vo as in Theorem [TJ To 
motivate our way of obtaining the Bessel equation as an approximating equation, 
consider the model operator 

P{x,hD) := -h 2 d 2 x + (x)- 2 

It is tempting to introduce the Bessel operator 

P (x,hD) := -h 2 d 2 x + x~ 2 
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which should be a good approximation for large x. The problem here is that even 
though the error decays like x~ 4 it is not small compared to H unless x > ftr* . Since 
we need to be able to send h and E to zero independently, such an approximation 
is useless for the case were E is small but fixed and K — > 0. To idea behind our 
reduction to the Bessel equation is essentially to let (x) be a new independent 
variable. The reader will easily see that this is precisely what Lemma @] does (in 
addition, we scale out E and fix the turning point to lie at 1). 

Lemma 12. For any < E < Eq the following holds: f(x) is a smooth solution of 

-h 2 f"(x) + V(x)f(x)=Ef(x) on x>E~i 

iff 0(f) = (j){^E) := n(£,E)^f(E-iy(Z,E)), with ft as in JX2| 7 is a smooth 
solution of 

(A.3) - fiV(0 + [r a (l - ^ 2 /4) - 1]0(O = h 2 W o (t,E)0(O on f > fi(-E) 
with a potential Wq satisfying 

(A.4) |0|9fWb(£, £)| < C M £- fe r 3 ~* 

/or > 0. 

Proof. Under the change of variables g(y) = f(E~?y) the following equations are 
equivalent, with V\(x) = — j(x)~ 2 : 

-h 2 f"{x) + (V (x) + h 2 Vx(x))f{x) = Ef(x) 

-h 2 g"{y) + (E^VoiE-iy) - l)g(y) = -h 2 E- 1 V 1 {E" 1 -y)g(y) 

Now let f = £(y,E) be as in Lemma 0] and set tpiO — 9(y(£iE))i or equivalcntly, 
ip(£(y,E)) — g(y). Then, with [i as in Lemma [TT1 

- »V'(0 + (d^,E)) 2 (d vv 0(y(^E),E)^(0] + (r 2 - 

(a.5) = -tfV'CO + W(0] + (r 2 - i)v(0 

= -h 2 {d iV {t E)fE- 1 V 1 {E-^, 
Let ft be as in Lemma [TT] In view of (|A.5[) . := ftiip satisfies the equation 

(a.6) -»V(o + (r 2 - iMO = » a w&£O0(O 

with 

The second part here involving ft decays like f -4 , whereas the first only decays like 
f ~ 2 . We need to extract this leading order decay: the asymptotic expansion 

M0 = -^ + O(C 3 ) as f^oo 

and Lemma HI imply that 

(dsy&E^E-^iE-iy&E)) = ~ + C 3 V r (t E) 

where 

\& E %V r {£,E)\ <C u E- k C l 
In view of (|A.6|1 . this yields equation (|A.3|) and we are done. □ 
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A fundamental system {^v^}? =1 of the homogeneous form of (|A.3j) . i.e., 

(a.?) - h 2 <p"{o + [r 2 (i - h 2 /A) - i\m = 

is given in terms of Hankel functions: 
or, equivalently, by the Bessel functions 
with Wronskian 

7T 7T 

Hence, the forward Green function of (|A.7|) is 

G n (tO ■= J[^(0C(O-C^')C(0]XK<f'] 

and thus a basis {<j>j,n}j=i °f (| A.3|) is given by the Volterra equation 

o l (e,f)^o(r,£)<MO^ / 

that one now needs to solve. This of course requires a thorough understanding of 
the behavior of Jn(n£) and Y n (n£) for large n on intervals of the form £ > £o > 
where < £o "C 1 is fixed, see [T] and [15]. We leave it to the interested reader to 
pursue this direction. 
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